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An Improved Limit on the Rate of the Decay 
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We report results of a search for the lepton-family number violating decay K + — > 7r + /i + e~ from 
data collected by experiment E865 in 1996 at the Alternating Gradient Synchroton of Brookhaven 
National Laboratory. We place an upper limit on the branching ratio at 3.9 x 10 -11 (90% C.L.). 
Together with results based on data collected in 1995 and an earlier experiment, E777, this result 
establishes a combined 90% confidence level upper limit on the branching ratio at 2.8 x 10 -11 . We 
also report a new upper limit on the branching ratio for tt° —> /j, + e~ of 3.8 x 10~ 10 (90% C.L.). 

PACS numbers: 13.20Eb, 14.40Aq, 11.30RV 



Despite the success of the Standard Model in describ- 
ing elementary particle physics, several key issues remain 
unresolved. One is the gauge-hierarchy problem, i.e., 
that the scalar Higgs field, introduced to give mass to 
the W and Z vector bosons, suggests a particle of mass 
f=a 100 GeV/c 2 while its renormalization implies a mass 
scale at least ten orders of magnitude larger. Theoreti- 
cal extensions to the Standard Model, such as Horizontal 
Gauge models Technicolor g], and Supersymme- 

try Q , were developed primarily to address the gauge hi- 
erarchy problem. These extensions permit lepton-family 
number (LFN) non-conserving decays. Observation of 
such phenomena would indicate new physics: physics be- 
yond the Standard Model. 

To test these theories in the kaon sector several ex- 
periments have recently been performed: K\ — > /i^e^ 
@, K° L 7r°/i ± e T |]'and K+ -» tt+ fi+e' {K n(us ). In 
the mid and late 1980s an experiment at the Brookhaven 
AGS, BNL E777, searched for the K vlie decay. Prior to 
that experiment the upper limit on the branching ratio 
for Knpe w as measured at CERN to be 4 x 10" 9 @ . E777 
reduced that limit to 2.1 x 10~ 10 (90% C.L.) @. In 1992 
a more sensitive experiment, BNL E865, continued the 
search for the K ntie mode. In the 1995 running period 
E865 achieved the same limit as E777 (10. In this pa- 
per, we report new results, from data collected in 1996. 

The E865 detector system and trajectories from a sim- 
ulated i^Tr^e decay are shown in Fig. |l]. Some of the de- 
tails of the detector have been previously discussed in the 



context of an associated measurement of the branching 
ratio K+ -> vr+e+e" fl|. 
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FIG. 1. Plan view of the E865 apparatus. 

The apparatus resided in an unseparated beam directly 
downstream from a 5 m long vacuum decay volume. With 
10 13 protons impinging on a 10 cm long Cu production 
target, a secondary unseparated 6 GeV/c beam was pro- 
duced containing about 10 8 K + and 2 x 10 9 ir + and pro- 
tons per 1.6 second AGS pulse. All detector elements 
were either desensitized or removed from the beam re- 
gion. A dipole magnet at the exit of the decay region 
approximately separated K + decay products by charge 
(negative to beam left, positive to beam right) and re- 
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duced the charged particle background originating up- 
stream of the detector. Proportional chamber packages 
(P1-P4), each containing four planes of chambers, were 
arrayed on either side of a second dipole magnet to form a 
momentum analyzing spectrometer system. The momen- 
tum resolution of this configuration was up as 0.003P 2 
GeV/c, where the momentum of the decay products, P, 
ranged from 0.6 to 4 GeV/c. 

Correct particle identification (PID) was of critical im- 
portance in reducing backgrounds. The first elements 
of the PID system were atmospheric pressure Cerenkov 
counters upstream and downstream of the spectrometer 
magnet. The left sides of each, C1L and C2L, were filled 
with hydrogen gas to detect e~ and not n~ , and had a 
light yield of as 1.7 photoelectrons (p.e.) for e~ . To re- 
ject e + , the right sides (C1R and C2R) were filled with 
methane with a light yield of about 4.5 p.e. for e + . The 
two sides were separated by a thin membrane. In order 
to reduce beam gas interactions, closed tubes of hydrogen 
gas were placed in the beam region of C1R and C2R. 

The second PID element was a Shashlyk-style electro- 
magnetic calorimeter (EM cal) [ p^| . This device con- 
sisted of 600 modules, each 11.4 by 11.4 cm 2 in cross 
section and 30 cm (15 radiation lengths) along the beam, 
arrayed 30 horizontally by 20 vertically with 18 modules 
removed in the beam region. The approximate resolution 
of the array for electrons was ue/E = 0.09 /y/E, where E 
is the electron energy measured in GeV. Typical energy 
deposition of a minimum ionizing particle was 250 MeV. 

The third clement was a muon range stack with 24 
planes of proportional tubes, with tubes alternately ori- 
ented horizontally and vertically, and plates of steel 
placed between each pair of planes. The steel thickness 
was 5 cm between the first eight pairs of planes and 10 
cm between the last four. 

Trigger hodoscopes were located directly downstream 
of PI (D hod), immediately upstream of the calorime- 
ter (A hod), and between the eighth and ninth pairs of 
proportional tubes in the range stack (B hod). 

The first-level trigger (TO) selected three charged par- 
ticle tracks based on hodoscope and calorimeter hit pat- 
terns consistent with kaon three body decays, and with 
at least one particle on each side of the apparatus. The 
w/j,e trigger added a signal from B hod, and at least 0.25 
p.e. from C1L and C2L. The TO trigger, downscaled by 
10 4 , also served as the trigger for our normalizing pro- 
cess, K + — > 7r + 7r + 7r~ (K T ). Satisfaction of any requested 
trigger, roughly 700 times per machine pulse, resulted in 
all information in the various data buffers being read out 
to a Fastbus based data acquisition system in about 100 
/is. 

This paragraph describes the most important offline se- 
lection criteria ]l3| . Common requirements for the K v/ie 
and the K T normalizing mode included a vertex formed 
from three reconstructed tracks, one negatively and two 
positively charged, with total vector momentum consis- 



tent with the beam phase space. Events with an extra 
(accidental) track on the left, not on the vertex, capable 
of making a C1L or C2L signal were removed. For K^^e 
candidates, additional cuts were required. Background 
from 7T° — ► 7e + e~ (Dalitz) was suppressed by removing 
events with M ee (effective mass of the e~ with cither 
positive particle interpreted as e + ) < 50 Mev/c 2 . Elec- 
tron PID required a signal in both C1L and C2L with 
corrected timing within ± 4 ns, and energy deposited in 
the calorimeter divided by the particle measured momen- 
tum (E/P) to be 1.0 ± 0.2. Pions were required to have 
C1R and C2R < 1.2 p.e., E/P < 0.85, and to not have 
penetrated the range stack to a depth of their expected 
ionization loss range. Muons were required to have C1R 
and C2R < 1.2 p.e., less than 450 MeV deposited in the 
calorimeter, and range stack penetration depth consis- 
tent with their range. High rates on the right, mostly 
at the one p.e. level, from high momentum muons and 
charged particle scintillation, prevented use of lower ir 
and /i thresholds in C1R and C2R. Table I summarizes 
the final PID efficiencies and probabilities of misidentifi- 
cation ]l3[ . These measurements were made with parti- 
cles of known identities from K T and from K^i and 
with a subsequent Dalitz decay of the 7r°. 
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0.873 ±0.002 








0.546 ±0.003 



TABLE I. PID efficiencies and probabilities of 
misidentification. The symbol "— denotes "identified 
as". 

A likelihood analysis was used to evaluate the proba- 
bility that selected events fit particular hypotheses, e.g., 
K^^e, K Tl accidental. In this method distributions de- 
rived from data were used as probability density func- 
tions (PDFs). These distributions included vertex and 
track quality, reconstructed beam parameters, PID de- 
tector responses and timing, and the invariant mass of 
the three decay products, to mention the most impor- 
tant. The extremes of these distributions were cut to 
allow the survival of about 95% of the events for the re- 
spective hypotheses. The resulting PDFs were then the 
templates to determine the probabilities that the vari- 
ables in a given event originated from the hypothesized 
mode, say the i th mode. The logarithms of these prob- 
abilities were added to form the joint log-likelihood (Lj) 
for the i th decay hypothesis. In the case of K ntie the 
PDF for the invariant mass of the decay products was 
generated by Monte Carlo simulation. All other PDFs 
were generated from data. 

An example of the final L T distribution is seen in Fig. |^ 
where we display data and Monte Carlo simulated K T 
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events. The 10% and 20% points on these plots represent 
likelihood values for which the probabilities of finding 
smaller likelihood are 10% and 20%, respectively. 



ane" response in C1L. The probability of such misiden- 
tification could then be calculated, and the number of 
such events normalized to the measured number of K T 
and Dalitz events in the total data sample. 
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FIG. 2. K T log-likelihood comparison between data and 
Monte-Carlo events. The vertical lines (solid for data and 
dash for MC) show the 10% and 20% likelihood points. 

Candidate K v(le events were first subjected to the cuts 
described above, and for those that survived the was 
determined. The results of this analysis are presented in 
Fig. as a scatter plot of h nfie vs. invariant mass of the 
decay products. The top plot is data where the invariant 
mass cut was increased from its nominal value of 3cr to 6<r 
for display purposes. The bottom plot shows simulated 
Kiltie events, where a vector interaction for the decay 
was assumed in the simulation. Only seven data events 
have survived the cuts with a L T/je greater than -170. 
Three of these are within the 3a accepted M nlie region, 
i.e., between the horizontal lines. The three data events 
which pass all cuts have probabilities of 13% (79%), 2% 
(52%), and 0.5% (64%), respectively, of being consistent 
with a K^fie (accidental) hypothesis. 

The three most probable sources of background events 
are K r , Dalitz, and accidentals. Contributions to back- 
grounds from the K T and Dalitz modes were from 
misidentification of particle species and incorrect track 
reconstruction, while accidental backgrounds were from 
events with correct PID, but for which the three parti- 
cles did not originate from a single decay and occurred 
accidentally in time. 

Estimation of the number and L w/je distribution for K T 
and Dalitz events was accomplished using events with 
correct PID, and replacing their detector response values 
with those corresponding to K vlie . This replacement was 
made by selecting randomly from a library of measured 
responses for events where the respective detectors gave 
incorrect PID, e.g., the C1L response of a ir~ from a K T 
event was replaced with one for which a known 7r~ had 
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FIG. 3. Scatter plot of E865 1996 K^e data (top) and 
Monte Carlo (bottom). The abscissa is the log-likelihood of 
the reconstructed events under the K vtie hypothesis, the or- 
dinate is the invariant mass of the detected particles. The 
horizontal lines demark the 3<r mass region. 

The number and h nfie distribution of accidental events 
was estimated using events for which the variable de- 
scribing the rms deviation from the mean time of all 
participating counters (T rms ) was more than 3 standard 
deviations, but which otherwise satisfied all K^^e cuts. 
Estimation of the relative probability that such events 
would have an acceptable T rms was made by evaluating 
the T rms distribution for events with total momentum 
greater than 6.5 GeV/c, i.e., events which are primarily 
accidental. The T rms distribution for high momentum 
events also gave the T rms PDF for accidental events with 
acceptable values of T rms . The latter was used to ran- 
domly replace the T rms value for accidental events de- 
scribed above in forming the full L^ Me distribution for 
accidentals. 

The Ljr^e distributions for K^^e, Dalitz, and accidental 
modes are shown in Fig. ^ [pj]. The estimated number 
of such events that would pass all K ntle cuts and appear 
in Fig. U within the accepted M x(ie region is 2.6 ± 1.0: 
0.06 ± 0.03 K T events, 0.1 ± 0.1 Dalitz events, and 2.4 ± 
1.0 accidental events, in good agreement with the three 
events observed. 

With the distributions for K-^^, and those of 

the most dominant background modes weighted in rel- 
ative proportion as determined above, a \ 2 function for 
Poisson-distributed data was minimized to determine the 
most probable number of K^^ e and total background 
events in the data distribution of Fig. ||. Those num- 
bers were 0.0 and 3.0, respectively, consistent with our 
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estimation of 2.6 background events. To determine the 
90% confidence interval for this result, the Frequentist 
approach was used Jl5| , with the number of background 
events assigned to be 2.6. Including the uncertainty in 
the assigned background level, the result of that analysis 
was that the expected number of K^^e events in our data 
sample is less than 2.5 at the 90% confidence level. 
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FIG. 4. Kwpe likelihood distributions for signal, Dalitz, and 
accidentals. The histograms are smoothed to compensate for 
the lack of statistics. The vertical scales for Dalitz and acci- 
dental distributions are normalized to their respective number 
of expected events. 

The upper limit on the K^^e branching ratio, normal- 
ized to the K T branching ratio, is calculated according to 
the formula: 



BR(irf±e) < BR{ititit) 



N (it fie) A^n 
N(nim) A nfie 



C, 



where BR denotes the noted branching ratio, BR(mnr) = 
0.0559 ± 0.0005 AT^e) = 2.5, the 90% C.L. num- 
ber of signal events; N(iririr) = 2.19 x 10 10 , the number 
of K T events adjusted for downscale factors; A repre- 
sents the geometrical acceptance of the detector system 
for the specific decay mode, A^utt IA vfxe = 1.64 ± 0.02, 
with 0.01 contribution from systematic uncertainty; and 
C = 3.78 ± 0.08, the product of correction factors ac- 
counting for efficiency differences between the two modes. 
The bulk of C, 3.15, is the reciprocal of the product of 
the ir + ,(i + ,e~ PID efficiencies shown in Table I, while 
the remaining 1.20 results from acceptance differences 
between K nfJie and K T due to cuts E^j. 

Employing these factors, we set a limit on the branch- 
ing ratio BR(wfie) < 3.9 x 10~ n (90% C.L.). Com- 
bining this result with data collected in 1995, BR < 
2.1 x 10" 10 JIB, and the E777 experiment, BR < 
2.0 x 10~ 10 yields a new upper limit of the branching 
ratio for K vlle of 2.8 x 10~ n (90% C.L.). 

This branching ratio implies that an intermediate bo- 
son in models described by a horizontal gauge interaction, 
e.g., Ref. with purely vector coupling and strength 



equal to that of the weak interaction, would have a mass 
greater than 60 TeV. 

Since the process 7r° — * /i + e~ would be observed in our 
data through K + — ► 7r + 7r°; tt° — > /i + e~, we also set an 
upper limit on its branching ratio. The only candidate 
events are the three discussed above, but their M Me val- 
ues of 0.226, 0.282, and 0.332 GeV/c 2 are too far from 
M^o for the events to have originated from tt° decays. 
We thus place an upper limit on the expected number of 
events at 2.44 (90% C.L.). We again normal- 



ize to the K. 

Attttti /A v 2(fj,e 



mode with the ratio of acceptances being 
= 4.07±0.02, the factor C=3.22±0.07, and 
the branching ratio for K V 2 is 21.16% jL6|. Combining 
these factors gives an upper limit on the decay branching 
ratio for vr° -> fi+e' of 3.8 x lO" 10 (90% C.L.), compared 
with the Review of Particle Physics limit of 1.72 x 10~ 8 
for the combined tt° — > {p. + e~ + fi~e + ) decays [ jillf . 
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